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1.  OBJECTIVE 


The  objective  of  the  propose  1  work  is  to  develop  a  computer  model  capable 
of  predicting  the  distribution  of  Bolid  particle  constituents  during  the 
processing  of  solid  propellants. 

2.  INTRODUCTION 

In  order  to  control  the  propellant  burning  rate  and,  hence,  the  propul¬ 
sive  performance  of  a  solid  rocket  motor,  it  is  essential  to  create  a  desired 
particle  distribution  throughout  a  solid  propellant.  It  is  therefore  desir¬ 
able  to  determine,  via  computation,  the  distribution  of  particulates  during 
casting  of  solid  rocket  fuel.  The  present  program,  fundamental  in  nature, 

>  -eared  toward  initiating,  understanding,  and  modeling  the  mechanisms  that 
control  the  distribution  of  solid  particle  constituents  during  the  processing 
of  a  solid  propellant. 

The  propellant  casting  can  be  considered  as  a  multiphase  flow  process 
with  a  dense  suspension  fluid  in  creeping  motion.  In  such  a  case,  the  shear 
resistance  (shear  lift  phenomenon),  relaxation  pher >menon  due  to  particle- 
particle  Interaction,  particle-fluid  interaction,  and  particle-wall  inter¬ 
action  predominate  even  though  the  fluid  phase  (binders)  may  be  serving  the 
main  function  of  transporting.  Even  if  we  start  with  a  uniform  mix  of 
particles  of  various  sizes,  non-uniformity  may  result  from: 

(1)  particle-fluid  interactions  alone  -  small  particles  due  to 
shear  lift  force  tend  to  move  away  from  the  wall  faster  than 
large  particles. 

(2)  particle-particle  interactions  alone  -  small  particles  tend  to 
act  as  a  fluid,  exerting  a  shear  lift  force  on  the  large  par¬ 
ticles.  The  large  particles  tend  to  move  away  from  the  wall 
faster  than  the  small  particles. 
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(3)  the  boundary  condition  at  the  wall  where  the  fluid  velocity 
must  be  zeto,  but  the  particle  velocity  may  have  a  finite 
valie  (slip  motion).  Thus,  particles  tend  to  migrate  toward 
the  center  of  the  passage  or  away  from  the  wall. 

These  expected  facts [1]  are  in  agreement  with  observation  by  Bradfield[2]  of 
WRE  (Weapons  Research  Establishment). 

In  treating  the  present  problem  of  transient  flow  with  solids  of  various 
particle  sizes  in  a  liquid  suspension,  we  need  to  consider: 

(1)  shear  lift  phenomenon, 

(2)  particle-particle  interaction, 

(3)  particle-fluid  interactions,  and 

(4)  particle-wall  interactions. 

The  shear  lift  phenomenon  of  solid  particles  in  a  fluid  was  formulated  by 
Saffman[3],  studied  experimentally  by  Segrfc  and  Silberberg[4] ,  and  correlated 
by  Soo[5],  including  the  resultant  concentration  in  steady  or  depositing 
flow[6].  The  effect  of  particle  diffusivity  was  found  to  be  important.  Segrd 
and  Silberberg[4]  particularly  identified  that  in  laminar  pipe  flow  of  a  sus¬ 
pension  of  neutral  buoyancy,  the  maximum  concentration  of  particles  tends  to 
occur  at  2/3  radius  from  the  pipe  axis,  that  is,  toward  the  wall,  than  toward 
the  center.  Many  of  these  basic  relations  concerning  the  distribution  of 
particulates  in  a  fluid  need  to  be  synthesized  and  formulated  for  the  present 
system. 

It  should  be  noted  that  very  little  work  has  been  done,  both  experimen¬ 
tally  and  theoretically,  v/a  dense  suspension  systems. 

3.  MATHEMATICAL  MODELING 

3.1  Multiple  Velocity  Field  Model  or  Multifluid  Model 

Use  of  a  multiple  velocity  field  model  calls  for  a  set  of  formulations 
including  continuity  and  momentum  equations  of  phases,  or  in  the  present  case. 
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components*  The  continuity  equations  are 


36  ,p  , 

— +  7  •  (VA)  -° 

36  jP  ^ 

iP+7  •  (<yA)  -o 

3e^p  ~ 

—  +  7  •  (*3P A)  ■  0 


(1) 


(2) 


(3) 


where  0j,  p  ^ ,  _Uj  refer  to  the  volume  fraction*  the  material  density*  and  the 
velocity  of  the  binder;  subscripts  2  refer  to  the  aluminum  powder  and  3  refers 
to  the  ammonium  perchlorate  powder.  For  the  example  of  6j  •  0.2352,  02  » 
0.1269,  and  63  ■  0.6379,  noting  that  0^  +  62  +  63  ■  1.  The  momentum  equations 
are 

aOjPjJLi 


ir-1  +  7  •  ( V ASi)  -  -  017P  +  7  •  li.  +  VA  -  ki2(Si  - 

-  k13  (Oj  -  a,) 

8e2p2^2 

-at - +  7  *  (02p2^2)  "  "  e27P  +  7  •  l2m  +  VA  '  *21^2  ’  *0 


(4) 


-  k23  (U2  -  U3) 


(5) 


863p3^3 


3t  +  7  ^03P3^3^3^  “  “  037P  +  7  *  53m  +  03p3^3  ~  K31^-3  ~ 

"  ^J2  CJL3  -  ^ 


(6) 


where  P  is  the  pressure,  t  is  the  shear  stress  of  component  1  in  the 
mixture,  Jjj  is  the  field  force  per  unit  mass  on  component  1  and  may  include 
that  due  to  the  shear  lift  effect*.  Kj2  is  the  Interfacial  momentum  transfer 
coefficient  Including  drag  between  phase  2  and  phase  1,  etc.,  where 


*The  magnitude  of  shear  lift  here  is  seen  not  to  be  influenced  by  the  rotation 
of  the  particles [3] . 


K12  "  9  IP  1F12»  etc 


Fj2  is  the  inverse  relaxation  time  for  momentum  transfer  from  phase  2  to  phase 
1,  etc.,  and 

9  lp  1F 12  "  e2P2F21  or  K12  "  K21 


etc.  from  action  and  reaction.  With  correct  boundary  conditions,  Eqs.  1  to  6 
are  solved  for  r-n  isothermal  system  to  determine  the  volume  fraction  and 
velocity  distribution  of  phases.  Transport  properties  are  needed  to  determine 

T  lm  etc*  an<*  f12  etc*  f12»  f21*  f13»  an^  F31  arise  from  fluid-particle 
interaction.  F23  or  F32  arises  from  particle-particle  interaction.  These 
quantities  depend  on  the  properties  of  materials  and  operating  conditions; 
non-Newtonian  behaviors  are  expected  for  the  present  system. 

For  given  initial  conditions,  pertinent  boundary  conditions  for  the  above 
equation  for  flow  through  a  pipe  of  radius  R  Include 
r  -  R,  Ulz  -  0 


3U2* 

U2z  "  "  L21  3r 


R 


au3 

U3z  "  -  L31  ~ 


for  the  axial  velocities,  where  L21  is  the  Interaction  length  of  particle  to 

/  2\1/2/  /  2\U2 

fluid,  leading  to  slip  motion  and  ■  ((AuJjj)  /  F21  w^ere  \Au2l/  *8 

the  relative  intensity  of  motion  of  phase  2  in  1,  and  F21  can  be  large  for 
small  particles  in  a  viscous  fluid.  A  limiting  case  will  be 

au2z/ar  -  au3z/ar  -  o. 

Since  the  particulate  material  finally  set  in  their  place  by  solidification 
rather  than  by  deposition,  the  boundary  condition  for  the  volume  fraction  of 
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particles  at  the  wall  la  given  by  16}  (k  -  2,  3) 

39k 

°k«9r  g'W'kl 

VbHa  Tja  la  defined  according  to  the  multiphase  formulation,  Its 
determination  for  the  present  system  la  complicated  because  of  a  dense 
suspension.  Based  on  the  theory  of  dense  suspensions [7 ] ,  is  expected  to 
be  greater  than  t j  (pure  binder).  In  general  xa  of  the  mixture  is  strongly 
influenced  by  the  perchlorate  powder  (63X  by  volume).  Unlike  the  case  of  a 
dilute  suspension,  t T2n*  ^3n  are  not  readily  determined  at  this  time.  (In 
the  case  of  dilute  suspensions  of  2  and  3,  Tj,  t2n  ■  D2m62p2»  T3m  " 

D3m®3P3»  °2n  the  dlffusivlty  of  particles  2  In  the  mixture.  Likewise,  D3m 
is  the  dlffuilvity  of  particles  3  in  the  mixture.)  In  the  present  case,  it 
suffices  to  say  that  T2a  =  W  2n^-2*  where  p2m  tlie  viscosity  of  phase  2  in  the 
mixture,  Aj  is  the  deformation  tensor  of  the  motion  of  phase  2,  and  Is 
related  to  D2a  according  to  the  relation  of  a  dense  suspension [7] , 

Computations  based  on  the  multlfluld  model  calls  for  simultaneous  solu¬ 
tion  of  Eqa.  1  to  6  with  pertinent  boundary  conditions  and  accurate  transport 
properties  F2j,  *31 »  *32  Ujn»  M2m»  an<i  p3m*  The  latter  are  not  readily 
computed  or  measured;  often  the  best  we  can  manage  is  the  viscosity  of  the 
mixture.  The  closeness,  though  different,  of  the  phase  velocities  also 
suggests  that  the  mixture  velocity  may  be  sufficiently  representative.  Thus, 
it  leads  us  to  use  the  relatively  simple  single  velocity  field  model  as 
described  in  the  following  section. 

3*2  Single  Velocity  Field  Model 

Summing  Eqs.  4,  5,  and  6  gives  the  momentum  equation  of  the  mixture  as 

“atS  +  v  •  (’AfU  +7  •  ,  1 3  W4--2.K4  -sj 


-  VP  +  V  • 


1+  +  <>J- 


Ta  *  while  Internal  action  and  reaction  for  moaantua  transfer  cancel  each 
other.  The  mixture  density  and  velocity  are  defined  by 

pa  "  l  Vk 
m  1,2,3  *  K 


pA  '  i.L A 


while  the  field  force  is  given  by 


Pgfg  “  I  Vk^k 
a-m  l23  k  K-r. 

for  k  ■  1,  2,  3  u  in  Section  3.1.  For  a  dense  mixture,  the  velocity 
difference  (U^  -  Uj  is  saall  and  the  third  term  on  the  left  side  of  Eq.  7  can 
be  neglected. 

Equations.  2  and  3  can  be  modified  by  considering  the  continuity  equation 
of  the  mixture  obtained  from  summing  Eqs.  1,  2,  and  3,  or 


+  V  •  pJJu  -  3  . 

The  continuity  equation  of  species  k  (k  "  2,  3)  can  be  rewritten  as 

-?r  +  7  •  -v  • 

Since  “  Jim)  "  Jfc>  the  general  flux  of  phase  k,  we  have 

+  7  "~7  ’4  • 


*pk8k 

.  -  -  -  +  U  •  V  p,  9.  -  -  p.  0.  V  •  U  -  V  •  J.  =  ~r~ 

O  t  m  K  K  it  K  U»  K  a  E 

m 

In  terms  of  mass  fraction  c^  -  0kpk^Po>  E<**  **  becomes 
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9p.ck 


Ft  +7  '  IvAl  -  ‘7  -4 


(13) 


Subtracting  the  product  of  c^  and  Eq.  9  from  Eq.  13,  and  rearranging,  we  gat 
3cu 


or 


».» r  +  (>A- 7  «k  -4 


^  °k 

p»rr-“7  •  Jk  • 


<U) 


(15) 


For  nearly  constant  pa,  the  coaponent  continuity  equation  takes  the  fora 
^  Plr®k  t  \ 

Tt  ■  ~n  +  4  ‘  7  (“A)  ■  - 7  •  4  • 


(16) 


Either  Eqs.  5  or  6,  neglecting  Inertial  forces,  pressure  gradient,  and 
assuming  equal  particle  velocities  and  binder  velocicy  approximately  equal  to 
the  mixture  velocity,  gives 


e-p-f. 

,  -  -  +  v 

^  Fkl 
8  ,j>  ,.f 


ft?'*) 


vrr  -  \m  7  (»A> 

kl 


(17) 


for  drift  by  field  forces  and  diffusion  by  concentration  gradient*  This  is 
because  shear  resistance  in  a  suspension  arises  from  resistance  to  transport 
of  momentum  by  diffusion.  The  kinematic  viscosity  and  dlffusivity  D^,  are 
related  by 


evPv 

ft 7  4  ■  -P  4. 7  4  *  -  V VW 

kl  kl 


(18) 


with  the  correspondence  ot  (v^/F^j)  V  ~  -D^.  This  correspondence  serves 
to  explain  the  relation  between  the  diffusion  model  and  the  multifluid 
aodex.  It  is  noted  that  the  diffusion  flux  is  usually  derived  in  a  different 
manner [8],  Equation  16  now  reduces  for  nearly  constant  pm  to 


(19) 


rr*  *  4  •  7  KV  ‘  •  (“k. 7  (9k»k)  -  4Vk/FkiI 

which  is  the  diffusion  equation. 

The  diffusion  equation  renders  the  continuity  equation  of  phase  k 
independent  of  its  momentum  equation.  However,  once  8^  is  determined,  can 
still  be  calculated  from  its  momentum  equation  if  it  is  needed.  Often  a 
knowledge  of  the  distribution  of  8^  is  sufficient.  Equations  7,  9,  and  19  (k 
■  2,  3)  can  be  solved  for  P,  8j,  82,  63  and  for  the  following  boundary 
conditions: 

r  -  R.  U.  -  0 


-  D 


3pk6k 


km  a  r 


f 

“  pk®k  F 


k_ 

kl 


for  given  initial  conditions. 

It  is  recognised  that  F2j,  F3P  F32,  Djy,,  and  still  have  to  be 
determined  experimentally  for  accurate  prediction. 


4.  EXPLORATORY  CALCULATIONS 


In  the  process  of  forming  a  solid  propellant  motor,  several  different 
particulate  materials  are  mixed  together  with  a  polymer  binder.  The  batch  Is 
mixed  until  the  mixture  becomes  uniform  and  homogeneous.  The  homogeneous 
mixture  Is  then  passed  through  a  network  of  pipes  and  ducts  to  a  mold.  As  the 
mixture  flows,  the  components  of  the  mixture  begin  to  displace  relative  to  one 
another  by  shear  motion.  This  gives  rise  to  non-uniform  propellant  properties 
in  the  mold  and  hence  the  final  cured  motor. 

The  ability  to  calculate  partial  component  separation  from  a  homogeneous 
mixture  is  a  crucial  feature  which  must  be  present  in  the  mathematical  model. 
In  order  to  investigate  and  demonstrate  that  the  proposed  mathematical  model 
is  capable  of  simulating  this  separation  phenomena,  some  exploratory  calcula- 
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tions  were  made  for  representative  situations.  The  numerical  results  pre¬ 
sented  here  were  obtained  with  the  COMMIX  code [9, 10],  which  was  modified  to 
carry  out  these  calculations. 

4.1  numerical  Model 

Flow  into  a  6.35E-3  m  (0.25  in)  gap  between  two  parallel  plates  ar.d  into 

an  annulus  is  considered.  A  two-velocity  field  model  is  used  to  describe  the 

flow  and  component  distributions.  (Component  #1  represents  the  polymer  binder 

and  component  #2  represents  the  particles.  The  governing  equations  are 

discretized  by  the  finite  volume  technique.  The  two-dimensional  computational 

domain  is  partitioned  into  8  equal  partitions  across  the  gap  and  20  along  the 

flow  direction.  The  overall  length  modeled  (0.0254  m)  was  long  enough  so  that 

the  flow  would  become  fully  developed.  At  the  entrance,  the  mixture  is 

assumed  to  be  homogeneous  and  have  a  uniform  velocity  of  0.01  m/s.  The  binder 

(component  #1)  is  assumed  to  stick  to  the  wall  (vj|w  ■  0) ,  while  the  particles 

/3v  2  \ 

(component  #2)  are  assumed  to  have  a  free  slip  boundary  condition  (  ^  )  «*  0. 

\  w  / 

A  semi-implicit  time-marching  algorithm  was  used  to  solve  the  system  of 
equations.  By  marching  in  time  until  all  quantities  (velocity  components  and 
volume  fraction)  converged  to  one  part  in  100,000,  a  steady-state  solution  was 
reached. 

4.2  Flow  Between  Parallel  Plates 

The  first  problem  considered  is  flow  between  two  parallel  plates 
separated  by  6.35E-3  m  (0.25  in).  A  homogeneous  50-50  mixture  (by  volume) 
enters  the  gap  uniformly  with  a  velocity  of  0.01  m/s.  Other  characteristics 
are  shown  in  Table  1.  By  the  time  the  flow  reaches  the  exit,  a  fully 
developed  situation  exists. 

Figure  1  shows  the  fully  developed  velocity  profiles  for  the  binder  (v^) 
and  the  particles  (V2).  The  differences  in  the  velocity  distribution  must  be 
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TABLE  1.  Flow  Between  Parallel  Plates 


Gap 

6.35E-3  m  (0.25  in) 

Inlet  velocity  vj  ■  V2 

0.01  m/s 

Inlet  0  j  (binder) 

0.5 

02  (particles) 

0.5 

Pi  “  w2 

200  Pa-s 

Pi  ■  P2 

1000  kg/m3 

K12 

1.0E8  Pa-s/m2 

Re 

3.17  5E-4 

entirely  due  to  the  different  boundary  conditions  imposed  on  the  2  components 
because  the  material  properties  used  are  the  same.  This  was  verified  by 
running  a  problem  where  both  components  had  no  slip  boundaries  and  this 
yielded  the  fully  developed  velocity  profile  v^  shewn  in  Fig.  1  for  both 
components.  For  reference,  the  wsll  known  parabolic  velocity  profile  is  also 
shown  and  labeled  vfl  in  Fig.  1.  Note  that  the  velocity  distributions  of  both 
the  binder  and  the  particles  are  flatter  than  the  reference  case.  The  par¬ 
ticle  velocity  is  consistently  higher  than  the  binder  velocity  throughout  the 
cross-section. 

Figure  2  shows  the  volume  fraction  distribution  of  binder  (0j)  and 
particles  (02) •  It  i8  readily  apparent  that  a  partial  separation  of 
components  has  occurred.  The  particles  have  a  relative  maximum  concentration 
in  the  center  of  the  gap  while  the  binder  is  more  concentrated  near  the  wall. 

It  is  worth  noting  that  this  separation  phenomena  is  due  entirely  to  the 
difference  in  velocity  boundary  conditions  at  the  wall  for  the  two  components. 
4.3  Annular  Flow 

The  second  problem  considered  is  two-dimensional  axisymmetric  flow 
between  two  concentric  cylindrical  surfaces  forming  an  annular  region.  The 
problem  characteristics  are  summarized  in  Table  2.  The  inner  radius  is 
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TABLE  2*  Plow  Into  Annulus 


Rin 

0.01905  m  (0.75  in) 

Rout 

0.0254  m  (1.0  in) 

Gap 

6.35E-3  m  (0.25  in) 

Inlet  velocity  vj  ■  V2 

0.01  m/s 

0  j  (binder) 

0.2242 

02  (particles) 

0.7758 

Wl  -  p2 

200  Pa-s 

p  j  (binder) 

920  kg/m3 

P2  (particles) 

1950  kg/m3 

k12 

1.0E8  Pa-s/m2 

0.01905  m  (0.75  in)  and  the  outer  radius  is  0.0254  (1.0  in).  This  results  in 
a  6.35E-3  m  (0.25  in)  annular  gap.  Here,  we  have  used  a  homogeneous  mixture 
with  a  higher  particle  volume  fraction  (.77)  and  a  material  density  difference 
between  the  binder  (920  kg/m3 )  and  particles  (1950  kg/m3). 

Figure  3  shows  the  fully  developed  'elocity  profiles.  The  velocity 
distribution  when  both  components  have  no  slip  boundary  conditions  is  labeled 
v^ .  This  solution  is  similar  to  the  parallel  plates  solution  except  the 
central  peak  velocity  occurs  nearer  the  inner  surface.  The  particle  velocity 
(V2)  is  consistently  higher  than  the  binder  velocity  (vp  throughout  the 
cross-section.  Due  to  the  higher  particle  concentration  and  material  density, 
the  velocity  profiles  are  flatter  and  maximum  velocities  less  than  the 
corresponding  results  for  the  parallel  plates. 

The  volume  fraction  distributions  are  shown  in  Fig.  4.  The  particle 
volume  fractions  (62)  show  the  highest  value  near  the  middle  and  lower  near 
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Fig.  4.  Volume  Fraction  Distribution  in  Fully  Developed  Annular  Flow 


the  walls.  The  Interfacial  drag  coefficient  (K12  ■  1.0E8  Pa-s/m2)  Is  the  same 
as  used  In  the  parallel  plates  case.  In  addition,  two  other  runs  were  made 
with  K  ■  1.0E9  and  K  ■  7.0E10.  The  results  are  shown  in  Fig.  4  where  for  K 
greater  than  1.0E10  the  separation  becomes  very  small. 

5.  DISCUSSION  OP  EXPERIMENTAL  MEASUREMENTS  FOR  TRANSPORT  PROPERTIES 

The  following  transport  properties  are  needed  for  the  system  with  a 
binder  and  2  different  sices  of  particles  under  consideration  so  that  the  set 
of  governing  equations  with  appropriate  initial  and  boundary  conditions  can  be 
solved. 

(1)  ym  nev'.ds  to  be  measured  vs  shear  rate  and  time  at  the  processing 
temperature.* 

(2)  and  D3m  to  be  measured  if  possible. 

(3)  f2  and  ^3  shear  liJt  force  and  gravity  force,  can  be  computed 
(please  see  Section  5.3). 

(4)  F21  and  F31  can  be  estimated  (please  see  Section  5.4). 

It  is  recognized  that  the  transport  properties  are  a  function  of  composi¬ 
tion  of  materials  and  time  histcxy.  He  recommend  that  the  time  history  and 
temperature  dependence  may  be  included  in  some  measurements  if  this  can  be 
conveniently  done. 

5.1  Viscosity  Measurements 

It  is  recommended  that  viscosity  be  measured  with  the  following 
compositions,  with  respect  to  time  and  shear  rate: 

( 1 )  Binder  only 

(2)  Binder  +  particle  2  (p2)  with  p2  at  1/3  nominal  value. 

Binder  +  p2  with  p2  at  2/3  nominal  value. 


*It  is  assumed  that  t  is  a  function  of  velocity  gradient  and  mixture 
viscosity— Non-Newtonian  fluid. 
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Binder  +  p2  with  p2  at  full  nominal  value* 

(3)  Binder  +  particle  3  (p3)  with  p3  at  1/3  nominal  value. 


Binder  +  p3  with  p3  at  2/3  nominal  value. 

Binder  +  p3  with  p3  at  full  nominal  value. 

(4)  Binder  +  p2  +  p3,  both  p2  and  p3  at  nominal  value. 

The  viscosity  u  of  a  fluid,  In  general,  exhibits  shear  rate  dependence 
and  can  be  characterised  as  follows 


'I 

j 


where  K  is  an  empirical  constant  and  1 3 W/3  r  |  is  the  shear  rate.  The  exponent 
n  accounts  for  various  rheological  behaviors  of  the  fluid.  n  ■  0  for  a 
Newtonian  fluid;  n  >  0  for  a  shear  thickening  fluid,  that  is,  a  viscosity 
which  increases  with  shear  rate;  the  reverse  is  the  case  of  a  <  0,  a  shear 
thinning  fluid.  Depending  on  the  shear  rates  and  the  particle  slse,  internal 
friction  or  collision,  a  given  fluid  may  exhibit  both  ranges  of  behavior [ 11] . 
For  the  case  of  a  shear  thickening  fluid  (n  ■  1) 

K  -  C  62  a?  o 
UP  P 


where  Cp  is  a  constant  (*  1  for  spheres  having  elastic  collision),  9p  is  the 
volume  fraction  of  particles  with  radius  a[ 1 ] • 

5.2  Diffualvity  Measurements 

It  is  recognized  that  diffusively  cannot  be  obtained  by  direct  measure¬ 
ment.  Therefore,  this  information  may  be  obtained  indirectly  through  densi¬ 
tometry  of  particle  distributions  and  photoelasticity  of  stress  measurement  of 
slices  from  solid  propellant  specimens.  We  recommend  that  these  measurements 
be  made  at  the  compositions  as  outlined  in  (1)  of  viscosity  measurements. 
Experimental  data  is  needed,  but  for  a  shear  thickening  fluid  (n  ■  I),  the 
particle  diffusivlty  is  given  by 
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k  -  2,  3 


where  C4  is  a  coefficient  of  the  order  of  10 [7]. 

3.3  Field  Forces  fj  and  J3 

These  include  gravity  g,  and  shear  lift  forces.  The  former  is  negligible 
for  the  present  case,  and  the  latter  has  been  given  analytically  for  spheres. 
It  will  be  determined  with  semi-empiricism  in  the  course  of  computation  by 
validation  from  particle  density  distributions  and  shear  stress  distributions 
in  the  solid  propellant  if  available.  It  is  given  for  a  spherical 
particle [3, 7]  by 

fk  ■  Cj 3<6.46>/4.](pm/pk)  [v  (“)  ||f|] 

where  W  is  the  velocity  outside  the  boundary  layer,  AW  ■  W  -  W^,  the  velocity 
difference  between  the  fluid  and  the  particle,  a  is  the  particle  radius,  and 
( 3W/3 r]  is  the  shear  rate  of  the  layer.  The  coefficient  Cj  accounts  for  non- 
spherlclty.  Lifting  of  larger  particles  by  shear  motion  of  small  particles 
will  be  dealt  with  in  a  later  report. 

5. A  Inverse  Relaxation  Time 

Inverse  relaxation  times  are  predictable  for  spherical  particles.  For 
the  present  application,  iteration  with  empirical  coefficients  and  validation 
by  final  particle  density  distributions  cin  yield  semi-empirical  modifications 
for  predictive  purposes.  For  particle  2,  F21  is  given  by 

7  5  e2u 

*l  F  2a2[pz  +  (pm/2)] 

where  Cp  is  an  empirical  coefficient,  82  is  the  volume  fraction  of  particles 
2,  p  is  the  fluid  viscosity,  and  a  is  the  radius  or  characteristic  dimension 
01  the  partide[l]. 
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If  all  the  ltena  in  Sections  5.1  end  5.2  can  be  furnished,  they  will  be 
Boat  helpful.  The  order  of  priority  should  be 

(1)  Viscosity  measurements. 

(2)  Dens it one try  of  particle  distribution  of  sliced  solid  apeciaena. 

(3)  Photoelasticity  of  sliced  solid  specimens. 

The  mininua  measurement  would  be  flow  through  a  tube  and  aeaaure  the 
flowrate  and  pressure  drop,  from  which  the  viscosity  and  non-Newtonian  para¬ 
meters  can  be  deduced.  We  can  proceed  with  results  from  Item  (1)  alone.  The 
effort  toward  achieving  a  realistic  computer  program  will  be  facilitated  or 
reduced  by  ha. mg  Items  (1)  and  (2),  and  more  so  If  Item  (3)  is  available. 


1  i few  f r:  Ktii  i Ki  0.  k- 


A  continuum  approach  to  the  modeling  of  a  dense  suspension  has  been  taken 
as  opposed  to  discrete  particle  tracking.  Within  the  continuum  approach,  two 
formulations  have  beer,  identified:  the  multlfluid  model,  and  the  diffusion 
model.  While  the  multifluid  model  Is  more  complete,  there  are  more  unknown 
coefficients  associated  with  the  model  and  these  coefficients  need  to  be 
determined.  The  single  velocity  diffusion  model  is  computationally  more 
economical  and  it  involves  relatively  fever  unknown  coefficients  than  the 
multifluid  model.  Correspondence  between  the  two  models  has  been  pointed  out. 

A  crucial  feature  needed  in  the  mathematical  model  is  the  ability  to 
predict  partial  component  separation  from  a  homogeneous  mixture.  The  only  way 
the  single  velocity  diffusion  model  can  predict  partial  separation  is  from  the 
shear  lift  field  force  term.  The  magnitude  cf  this  term  must  be  determined 
experimentally.  Due  to  the  current  lack  of  detailed  and  reliable  experimental 
data,  approximations  were  made  for  the  interaction  terms  in  the  multifluid 
model.  Even  in  the  absence  of  a  shear  lift  field  force,  the  nultlfluld  model 
can  predict  partial  component  separation  by  having  different  component 
velocity  boundary  conditions  at  the  wall. 


In  summary,  our  preliminary  investigation  of  flow  modeling  during  solid 
propellant  processing  has  yielded  the  following  conclusions: 

1.  Meaningful  prediction  of  concentration  distribution  of  components 
can  be  obtained  from  computations  with  a  minimal  acquisition  of 
transport  properties  at  the  initiating  phase.  These  preliminary 
predictions  can  be  used  to  guide  experiments  which  are  urgently 
needed  to  quantify  the  transport  properties  and  to  validate  the 
mathematical  models. 

2.  Currently,  the  miltifluid  model  and  the  single  velocity  diffusion 
model  are  viewed  as  being  complementary.  While  the  multifluid  model 
can  be  used  to  gain  insight  Into  the  underlying  individual  physical 
mechanisms,  the  diffusion  model  gives  global  phenomenological 
behavior  of  the  system.  The  interrelations  between  the  two  mathe¬ 
matical  models  have  been  clarified.  It  is  anticipated  that  at  a 
later  date,  depending  on  availability  of  the  needed  experimental 
data  and  understanding  of  physical  mechanisms,  we  shall  select  one 
of  the  two  models  as  the  reference  predictive  tool. 

3.  Validation  of  computed  concentrations  can  be  made  by  sections  of 
solidified  models  or  checked  by  burning  rates  for  uniformity  of 
pressure. 

4.  Exploratory  calculations  have  shown  the  ability  of  the  lmiltlfluid 
model  to  compute  partial  separation  of  components  by  boundary 
condition  differences.  More  parametric  study  can  be  done  to  give 
insight  into  the  sensitivity  of  the  various  empirical  coefficients. 

5.  The  problem  under  investigation  is  important,  but  very  difficult. 
The  preliminary  results  from  -he  present  study  lays  the  foundation 
for  the  future  work  and  it  appears  that  useful  results  can  be 
obtained. 
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Langley  Research  Center 
ATTN:  MS-185  Tech.  Lib. 

Hampton,  VA  23665 

National  Aeronautics  &  Space  Admin. 

Lewis  Research  Center 
ATTN:  Lib(D.  Morris) 

21000  Brookpark  Rd. 

Cleveland,  OH  44135 

National  Aeronautics  &  Space  Admin. 

Lewis  Research  Center 

ATTN:  MS-501-5,  D.  A.  Petrash 

21000  Brookpark  Rd. 

Cleveland,  OH  44135 

National  Aeronautics  &  Space  Admin. 

Lyndon  B.  Johnson  Space  Center 
ATTN:  JM2/Tech.  Lib. 

Houston,  TX  77058 

National  Aeronautics  &  Space  Admin. 
Scientific  Technical  Info.  Fac. 

ATTN:  Accessioning  Dept. 

P.0.  Box  8757 

Baltimore  Washington  Inti  Airport,  MD  21240 
Commander 

Naval  Air  Dev.  Ctr. 

ATTN:  Code  8131 
Warminster,  PA  18974 

Commander 

Naval  Air  Sys.  Comd. 

ATTN:  NAIR-00D4-Tech  Lib. 

Washington,  DC  20361 
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Commander 

Naval  Air  Sys  Comd 

ATTN:  NAIR-320G,  Mr.  Bertram  P.  Sobers 
Washington,  DC  20361 

Commander  Officer 
Naval  Intel  Spt  Ctr 
Information  Svc  Div 
ATTN:  Doc  Lib 
4301  Suitland  Rd. 

Washington,  DC  20390 

Commanding  Officer 

Naval  Ord  Sta-Indian  Head 

ATTN:  Tech  Lib,  Code  4243C,  Henrietta  Gross 

Indian  Head,  MD  20640 

Superintendent 

Naval  Postgraduate  Sch. 

ATTN:  Code  1424-Libs  Dir 
Monterey,  CA  93943 

Chief  Naval  Research 

ATTN:  Dr.  Richard  S.  Miller,  Code  432 
Arlington,  VA  2221.7 

Chief 

Naval  Research 

ATTN:  R.  Junker,  Code  412 

Arlington,  VA  22217 

Commanding  Officer 
Naval  Research  Lab 
ATTN:  Code  6100 
Washington,  DC  20375 

Director 

Naval  Arsearch  West  Pasadena 
ATTN:  R.  J.  Mareus 
1030  E  Green  St 
Pasadena,  CA  91106 

Commander 

Naval  Sea  Sys  Comd 
ATTN:  SEA-09B312,  Tech  Lib 
Natl  Ctr  BG  3 
Washington,  DC  20362 

Commander 

Naval  Sea  Sys  Comd 

ATTN:  Mr.  Elgin  Werback,  SEA-62Z31B 

Natl  Ctr  BG  3 

Washington,  DC  20362 
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Commander 

Naval  Surface  Wpns  Ctr 
ATTN:  Acquisitions,  Code  E431 
Dahlgren,  VA  22448 

Commander 

Naval  Surface  Wpns  Ctr 

ATTN:  Code  E432,  S.  Happel,  Room  1-321 

Silver  Spring,  MD  20910 

Commanding  Officer 
Naval  Underwater  Sys  Ctr 
ATTN:  Tech  Lib  021312 
Newport,  RI  02840 

Commander 

Naval  Weapon  Center 
ATTN:  Code  343 
China  Lake,  Ca  93555 

Director 

Navy  Strat  Sys  Proj  Ofc 
ATTN:  Tech  Lib  Br  Hd 
Washington,  DC  20376 

Commander 
Ogden  ALC  (MANPA) 

ATTN:  Mr.  Anthony  J.  Inver so 
BG  1941 

Hill  AFB,  UT  84056 
Commander 

Radford  Army  Ammo  Plant 
ATTN:  SMCRA-QA 
Radford,  VA  24141 

Rockwell  Int'l  Corp. 

Rocketdyne  Div 
ATTN:  Tech  Info  Ctr 
6633  Canoga  Ave 
Canoga  Park,  CA  91304 

Rohm  &  Haas  Co. 

ATTN:  Scty  Off,  (Dr.  H.  M.  Shuey) 

723-A  Arcadia  Circle  1 

Huntsville,  AL  35801 

Commander 
SAALC  (SFTT) 

ATTN:  W.  E.  Vandeventer 
Kelly  AFB,  TX  78241 


1 


SRI  Int'l  Document  Ctr 

ATTN:  Classified  Doc  Svc,  (Dr.  Clifford  D.  Bedford) 

333  Ravensvood  Ave 
Menlo  Park,  CA  94025 

Talley  Industries 

ATTN:  Eng.  Tech  Lib,  (Kim  St.  Clair)  1 

P.0.  Box  849 
Mesa,  AZ  85201 

Teledyne  Ryan  Aeronautical 

ATTN:  Tech  Info  Svcs  (W.  E.  Ebner)  1 

2701  Harbor  Drive 
San  Diego,  CA  92101 

Thlokol  Chem  Corp 
Wasatch  Dlv 

ATTN:  Tech  Lib  (J.  E.  Hansen)  2 

Brigham  City,  UT  84302 

Thlokol  Corp. 

ATTN:  Scty  Off,  (D.  J.  McDaniel)  1 

P.0.  Box  241 
Elkton,  MD  21921 

Thlokol  Corp. 

ATTN:  Tech  Lib  (H.  H.  Sellers)  1 

Huntsville,  AL  35807 

TRW  Inc. 

Electronics  &  Defense  Sector 

ATTN:  Tech  Inf  Ctr,  Doc  Svcs  for  S/1930  2 

One  Space  Park 

Redondo  Beach,  CA  90278 

TRW  Inc. 

Electronics  &  Defense  Sector 
ATTN:  Tech  Inf  Ctr,  Doc  Svcs 

for  R.  C.  Reeve,  San  Bernardino  1 

>ne  Space  Park 
Redondo  Beach,  CA  90278 

USDRE  (PCA) 

ATTN:  0USDRE&E  (R&AT/MST) ,  (Dr.  Robert  J.  Heaston)  1 

The  Pentagon,  Room  3D1089 
Washington,  DC  20301 

United  Technologies  Corp 
Chemical  Systems  Div 

ATTN:  Tech  Lib  1 

P.0.  Box  358 
Sunnyvale,  CA  94088 
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United  Technologies  Corp. 

Research  Center 

ATTN:  Acq  Lib  (M.  E.  Donelly)  1 

400  Main  Street 

Eaat  Hartford,  CT  06108 

Conaander 

White  Sands  Missile  Range 

ATTN:  Tech  Lib  1 

White  Sands  Missile  Range,  NM  88002 


AMSMI-RD,  Dr.  McCorkle 
Dr.  Rhoades 

AMSMI-RD- RE,  Dr.  Hartman 
Dr.  Bennett 

AMSMI-RD-PR,  Dr.  Stephens 
Dr.  Wharton 
AMSMI-RD- PR-T,  Dr.  Alley 
Mr.  Allen 

AMSMI-RD-PR-P,  Mr.  Schultz 
AMSMI-RD-PR-M,  Mr.  If shin 
AMSMI-RD-PR-E,  Mr.  Maykut 
AMSMI-RD-CS-R 
AMSMI-RD-CS-T 
AMSMI-GC-IP,  Mr.  Bush 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

5 

15 

1 

1 


*C.S.  OOVTUMSST  PHINTHW  QITICJi  1968-530-081/80110 


DIST-10 


